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Abstract 

The corrosion of materials is an undesirable and costly process affecting many areas of technol¬ 
ogy and everyday life. As such, considerable effort has gone into understanding and preventing it. 
Organic molecule based coatings can in certain circumstances act as effective corrosion inhibitors. 
Although they have been used to great effect for more than sixty years, how they function at the 
atomic-level is still a matter of debate. In this work, computer simulation approaches based on 
density functional theory are used to investigate benzotriazole (BTAH), one of the most widely 
used and studied corrosion inhibitors for copper. In particular, the structures formed by protonated 
and deprotonated BTAH molecules on Cu(lll) have been determined and linked to their inhibiting 
properties. It is found that hydrogen bonding, van der Waals interactions and steric repulsions all 
contribute in shaping how BTAH molecules adsorb, with flat-lying structures preferred at low cov¬ 
erage and upright configurations preferred at high coverage. The interaction of the dehydrogenated 
benzotriazole molecule (BTA) with the copper surface is instead dominated by strong chemisorption 
via the azole moiety with the aid of copper adatoms. Structures of dimers or chains are found to be 
the most stable structures at all coverages, in good agreement with scanning tunnelling microscopy 
results. Benzotriazole thus shows a complex phase behaviour in which van der Waals forces play an 
important role and which depends on coverage and on its protonation state and all of these factors 
feasibly contribute to its effectiveness as a corrosion inhibitor. 


1 Introduction 

The oxidation and corrosion of metals affects many areas of industry, technology, and everyday life. 
Indeed, it plays a critical role in the failure of metal parts in engineering, electrochemical and catalytic 
devices, as well as of metal construction parts such as piping and roofing, to name but a few. The 
atomistic details of the processes of metal corrosion are still unclear, however an effective solution 
to the problem of corrosion has been found in the use of organic molecules as corrosion inhibitors, 
which work by interacting with metal/oxide surfaces and forming a protective fihrf 12 . Widely-used 
organic corrosion inhibition systems, whose effectiveness has been empirically verified, are amines or 
zinc dithiophosphates on steel and benzotriazole (BTAlU) on copper. However, the chemical details 
of the inhibiting action and the structure of the protective layer these molecules form against metal 
surfaces are still unresolved. In this work, state-of-the-art computational methods are applied to the 
BTAH/Cu copper system, which is the most studied both experimentally 4 and computationally 5 10 . 
The aim is to identify the structures BTAH forms on the copper surface and obtain insights into how 
they affect its function as a corrosion inhibitor. 

Most experimental studies on organic corrosion inhibitors are macroscopic studies involving the 
immersion of a metal specimen in a solution containing corrosive agents and organic inhibitors. In¬ 
formation about inhibition can be gathered during sample immersion, monitoring the evolution of the 
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material using electrochemical or spectroscopic techniques (see e.g. Refs. 11 13 ) or after immersion, 
by measuring the mass change of the specimen or the amount of metal solvated in the solution (see 
e.g. RefP). These studies reveal insight into the effectiveness of different chemicals. An alternative 
approach is to examine well-defined copper single-crystal surfaces under ultra-high vacuum (UHV) 
using X-ray absorption spectroscopy (XAS) and scanning tunnelling microscopy (STM). These tech¬ 
niques, used mainly by the surface science community, reveal structural details of the adsorbed system. 
The inhibitor molecules can be deposited onto the surface either from gas phase or via an aqueous 
solution (which is subsequently evaporated). Discrepancies are seen for different experimental condi¬ 
tions, with zigzag structures of flatly-adsorbed molecules having been postulated 11 13 (for molecules 
deposited from aqueous solutions), as well as a structure of vertically-adsorbed BTAs 15 (for evapo¬ 
rated molecules). In this work we only compare our results with experiments performed on evaporated 
BTAH molecules on copper substrates in UHV, since the computational work presented here consists 
of molecules adsorbed on a Cu(lll) surface in vacuum conditions. Previous studies have shown that 
BTAH deprotonates when adsorbed on copper surfaces, revealing that it can deprotonate not only by 
interacting with high-pH environments (its pK a is 8.4 at 25° C) but also through its interaction with 
the surface. BTA was seen to form nearly-upright organometallic surface complexes involving bonds 
between the azole nitrogen atoms and copper adatoms 15 18 . In particular STM studies showed the 
formation of strings of stacked dimers at low coverage and of more complex structures, also composed 
of BTA-Cu a d-BTA dimers, at high coverage. 

A number of computational studies, using density functional theory (DFT), have looked at the 
stability of BTAH and BTA structures on copper. Benzotriazole is a challenging molecule to study 
with DFT as it combines a strongly electronegative azole moiety, which preferably interacts with the 
surface through chemisorption, with a benzene-like ring which can interact with the surface via van 
der Waals (vdW) forces. Traditionally, vdW dispersion forces have been a challenge for DFT methods 
and are not accounted for in the most widely used exchange-correlation functionals, such as Perdew- 
Burke-Ernzerhof (PBE) 1 ^, that have been used in most previous studies of BTAH. Thus the role of 
vdW forces in this type of system is still largely unexplored. Isolated BTAH was found, in calculations 
using PBE, to chemisorb weakly to a Cu(lll) surface in an upright geometry, forming two N-Cu bonds 
via the triazole moiety 56 . Physisorption, with the molecule lying flat onto the surface, was found to 
yield a very small binding energy (~ 0.1 eV) with PBE 6 . The addition of an empirical van der Waals 
correction lead instead to a much stronger bond ( 0.72 eV)^. At higher coverages^, BTAH was 

found to form hydrogen-bonded (HB) chains with the molecules lying parallel to the surface. The 
dehydrogenated BTA molecule was found to bind strongly to the surface in an upright configuration 5 
when isolated, and to form Cu-bonded organometallic complexes of vertical or tilted molecules at 
higher coverages. However, the structure of the organometallic BTA-Cu ac j complexes is still debated. 
Chen and Hakkinen 10 found that deprotonated BTA-Cu ac j-BTA dimers are more stable than [BTA- 
Cu a d]n chains thus agreeing with the experimental results of Grillo et al . ^ and in disagreement with 
the DFT results of Kokalj et al. 8 . 

Here we report the most extensive DFT study performed to date and we explore in depth the 
importance of using vdW approaches in the study of adsorbate-surface interactions. Intramolecular 
interactions and coverage effects are also examined in detail; they are found to be particularly relevant 
for the systems formed by fully protonated BTAHs, with hydrogen bonding dominating in the low- 
coverage structures and vdW and electrostatic forces at high coverage. Moreover, the energetics of the 
dissociation process and of the formation of complex structures with Cu adatoms were investigated 
and linked to experimental conditions. The obtained structures and adsorption energy trends are 
comparable to the experimental results in UHV and a link with the effectiveness of the molecule as a 
corrosion inhibitor is discussed. 

The remainder of the paper is organised as follows. The computational methodology and set-up is 
presented in the next section (Sec. [2]), followed by the results for the protonated (Sec. [3]) and deproto- 
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nated (Sec. [4]) molecules adsorbed on Cu(lll). Finally, a discussion and conclusions are presented in 
Sec. 0 

2 Methodology and Computational Setup 

Calculations of inhibitor molecules adsorbed on copper surfaces were performed by means of DFT 
using the VASP code 20 23 . In the present work most results have been obtained with the optB86b- 
vdW functional, a modified version of the non-local vdW density functional 24 which explicitly accounts 
for dispersion based on the electron density. Indeed, optB86b-vdW has been shown to perform best in 
comparison to experiment for several adsorption problem s^ 25 ! 26 !, including the (relevant for this work) 
adsorption of benzene on Cu(lll) 27 . A number of other functionals were also considered for testing 
and comparison with previous results: PBE and a number of vdW-inclusive functionals (vdW-DF®!, 
optB86b-vdw! 28 } PBE-D2®! and PBE-TS®). For all functionals the calculated value for the lattice 
constant a and of the bulk modulus B are within 10% of the experimental value (see Table [l]), and in 
good agreement with previous theoretical results 28 . 



a [A] 

B [GPa] 

PBE 

3.64 

137 

vdW-DF 

3.71 

113 

optB86b-vdW 

3.60 

148 

PBE-D2 

3.57 

148 

PBE-TS 

3.55 

168 

Expt. 

3.49 

137 


Table 1 Lattice constants and bulk moduli for all the exchange-correlation functionals used in this work. All 
the values are in good agreement with the experimental values corrected for the zero-point anharmonic 
expansion® 1 1 and with previous computational results 2 -. 


Slabs of 3-6 layers were tested with optB86b-vdW (the functional chosen for the calculations 
of the BTAH/Cu systems), and the surface energy was seen to converge at a 4-layer thickness to 
11 _VdW = 0-102 eV/A 2 , in good agreement with the experimental value®! ( 7 ^( 111 ) = 0.112 
eV/A 2 ). In adsorption calculations the atoms in the bottom layer of the slab were kept fixed at the 
bulk positions and a vacuum of ^ 20 A in the non-periodic direction was added to prevent interactions 
between molecules in neighbouring cells. All calculations were performed using the PAW method 32 33 
and a kinetic energy cutoff on the plane wave basis set of 400 eV. The Brillouin-zone integration was 
performed using Monkhorst-Pack grids of (12x12x12) k-points for the conventional unit cell of copper 
bulk and of (12 x 12 x 1) k-points for a 1 x 1 Cu(lll) surface with the vacuum in the z direction. 
A dipole correction as implemented in the VASP code 34 35 was added in all cases, and calculations 
of the dissociated molecule (into BTA and H) in the gas-phase were spin-polarized. The inhibitor 
(BTAH and BTA) molecules were adsorbed onto the surface at a variety of coverages, between 1/16 to 
1/4 of a monolayer (ML), where 1 ML corresponds to 1 adsorbed molecule per surface copper atom. 
Experimental studies looking at the sub-ML regime were performed at 1/16 ML 18 and at 1/3 Mli®, 
close to the coverages studied in this work. A variety of hexagonal and orthorhombic cells were used 
and a large number of initial configurations have been chosen and optimised, in order to explore phase 
space as comprehensively as possible. 

The adsorption energies, E^ S: were calculated as follows: 


E^ds ( E S y S Egiab N mo \ X -E/ mo i ) /V mo i, 


( 1 ) 


where E sys is the total energy of the adsorbed complex, E s i a b is the total energy of the substrate (a 
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Cu(lll) slab with or without adatoms), E mo i is the total energy of the isolated molecule and N mo \ is 
the number of molecules in each cell. A negative value of indicates a stable adsorbed structure. 

The formation of complex systems, such as deprotonated BTAH on a copper surface with de¬ 
fects, includes endothermic (bond breaking and defect formation) and exothermic (adsorption of the 
molecules on the surface) processes. In order to account for all these energy contributions, an energy of 
formation, Ef orm , is defined here and used as an indicator of the likelihood of formation of the system. 
The following two stoichiometric equations are relevant to the formation of a system of adsorbed BTA 
molecules. In the first one, n BTAH molecules dissociate and adsorb on a Cu(lll) surface where m 
copper adatoms Cu a d are already present, to form the nBTA/Cu(lll)/mCu a d organometallic complex: 

nBTAH + Cu(lll)/mCu a d —» nBTA/Cu(lll)/mCu a d + H sys . (2) 

The second process is similar, however the copper surface is assumed to initially be atomically flat 
Cu(lll) and m copper adatoms Cu a d are formed from the bulk: 


nBTAH + Cu(lll) + mCubuik —^ nBTA/Cu(lll)/mCu a d + H sys . (3) 


In Eqs.[2p| the term H sys indicates either of two possible scenarios for the dissociated hydrogen atoms: 


they can either adsorb onto the copper surface or combine to form gas-phase hydrogen molecules. The 
formation energy, Ef orm , is defined as the difference in total energy between the final (right-hand side 
of Eqs. 2][3) and the initial system (left-hand-side of Eqs. [2][3]). Note that the adsorption energy E a d s 
of BTAH on clean Cu(lll) can be considered as a formation energy obtained from the stochiometric 
equation: 

nBTAH + Cu(lll) -A nBTAH/Cu(lll). (4) 


For certain overlayer structures, energy decompositions have been performed in order to estimate 
the strength of hydrogen bonds, of non-HB intermolecular interactions and of adsorbate-substrate 
bonding. These are performed by taking the optimised system and performing single point calculations 
of the individual system components. The strength of a hydrogen bond, Ehb, is obtained by performing 
single point calculations for isolated (gas phase) HB chains or dimers and for each gas phase BTAH 
molecule: 

E H b = (E^btah _ T e btah)/N, (5) 

N 

where E^ BTAH is the total energy in gas phase of the HB system of NxBTAH molecules, and Ebtah 
is the total energy of each BTAH molecule in the chain. Similarly, non-HB interactions are estimated 
using: 

Enon-HB = (EjJfxBTAH - 22 E B TAH ~ E H b)/N, (6) 

N 

where EJ^btah is the total energy obtained by removing the copper surface from the unit cell and 
performing a single point calculation, Ebtah is the total energy of each BTAH molecule in the system 
and Ehb is the strength of the HB interactions (if any) calculated using Eq. [5j The interaction energy 
of the molecules with the substrate is obtained using: 


Emol/Cu(lll) — (E S ys — E^btah — ^Cu(m))/N, (7) 

where E sys is the total energy of the adsorbed system, EJJfxBTAH and Ec u ( 111 ) the total energies obtained 
by performing a single point calculation after removing, respectively, the surface and the molecules 
from the unit cell. Energy decomposition schemes entail a certain level of arbitrariness®!, however they 
are also useful in providing a semi-quantitative understanding of the relative importance of hydrogen 
bonding versus adsorbate-substrate bonding. 
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The nature of the bonding between the adsorbates and substrate was investigated by looking at 
electron density difference plots. They were constructed by subtracting from the charge density of 
an optimised adsorbed system the charge density of the adsorbates and substrate. As for the energy 
decompositions, single point calculations are performed on the components i.e. they are kept to the 
structure they assume in the adsorption system. 


3 BTAH on Cu(lll) 


The adsorption of benzotriazole on Cu(lll) was studied first for a variety of coverages. The Cu(lll) 
surface was chosen since it has the lowest surface energy among the copper surfaces, and therefore it 
would be the most prevalent surface even on a copper polycrystalline sample. Moreover, it is the most 
widely studied surface experimentally and computationally. 

Although experiments in vacuum have shown that BTAH in contact with a copper surface depro- 
tonates^KS it is still useful to look at the adsorption of the fully protonated molecules for at least two 
reasons. First, the drive to deprotonation can depend on the environment and in acidic environments 
the molecule can be expected to be protonated. Indeed, a strong dependence on pH has been found 
on the effectiveness of benzotriazole as a corrosion inhibitor 37 (with worse inhibition in acidic media), 
suggesting that the molecular structure changes with pH and, thus, with the BTAH/BTA ratio present 
in the system. Second, the BTAH/Cu(lll) system is well-studied and it provides the possibility of 
validating the results obtained here with previous work. 

A large number of DFT calculations (summarized in Fig. ro was performed for BTAH on Cu(lll) 
considering molecules at different adsorption sites, different orientations and different coverages. In 


this section we are going to look first at the behaviour of isolated molecules (Sec. 3.1), and subsequently 
at supramolecular complexes (Sec. |3.2| ). In particular, we are going to highlight the role of hydrogen 
bonding, van der Waals dispersion and chemical bonding in the coverage dependence of BTAH on 
Cu(lll). 


3.1 Isolated BTAH on Cu(lll) 

The most stable adsorption configurations for a single BTAH on Cu(lll) at low coverage (1/16 ML) 
were established using both PBE and the vdW functionals mentiond in Sec. [2| Four stable structures 
were obtained and they are shown in Fig. [TJ Two upright configurations were found: the c Up’ structure, 
which binds to the surface via the N2 triazole atom and the ‘Side’ structure which binds via the N2 
and N2 triazole atoms. A physisorbed flat configuration (‘Flat’ in Fig. |T]) was also seen. Finally, a 
not-previously seen tilted configuration (Tilted’ in Fig. [T]) which combines chemisorption through two 
N-Cu bonds (as shown in Fig. [2]), with physisorption through 7r-bonding via the benzene-like ring, was 
identified. 

The comparison of their adsorption energies obtained with different functionals highlights the im¬ 
portance of choosing the most suitable methodology for the system under study. Indeed, the graph in 
^g. a shows that the relative adsorption energies of the four structures changes when they are opti¬ 
mised with different functionals. For example, the least stable structure with the PBE functional is the 
‘Flat’ configuration, since PBE is unable to describe physisorption. This is however one of the most 
stable structures for all the other functionals. The adsorption energy of the flat structure with PBE 
is E ac i s — —0.10 eV/mol, consistent with Ref. 38 . It is also in line with PBE calculations of benzene 
on Cu(lll) (E^ zene = —0.08 eV/mol) 28 which also showed almost no interaction of the molecule 
with the surface. When dispersion interactions are accounted for, the stability of the system increases. 
In addition, we find that the adsorption energy of the ‘Flat’ structure is linearly related to that of 
the benzene ring (Fig. [lj, and it was found to be between = —0.66 to E F ^ E_TS = —1.29. 

The nature of the interaction of the molecule with the surface, i.e. physisorption, was confirmed by 
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Figure 1 Top: adsorption energies for the Flat, Tilted, Side and Up structures calculated with five different 
functionals. They are plotted as a function of the adsorption energy of benzene on Cu(lll), calcu lated with 
the corresponding functional. The adsorption energies of benzene on Cu(lll) are taken from Refs.^EHl. 
Bottom: structures of the four lowest energy configurations for a single BTAH adsorbed on Cu(lll). Selected 
bond distances are given in A. The molecule and labelling scheme are shown on the left. 


looking at electron density difference plots, shown, for selected functionals, in Fig. [2] No interaction 
between the surface and the molecule is observed via the triazole ring, and only a small amount of 
charge rearrangement is observed for functionals which yield a stronger binding of benzene to Cu(lll) 
(PBE-TS and optB86b-vdW). The molecule-surface distance, when optimizing with functional which 
underbind the benzene ring (PBE and vdW-DF), is constant across the molecules (3.55 — 3.65 A). 
When functionals which strongly bind the benzene ring are used, a small amount of tilting is observed, 
in good agreement with Ref. 39 . 

The most stable configurations with PBE are the upright chemisorbed configurations (‘Up’ and 
‘Side’ in Fig. [TJ. Their PBE adsorption energies ( E ac js = —0.40 eV/mol for ‘Up’ and E^ds — —0.36 
eV/mol for ‘Side’) were found to be in excellent agreement with Ref. 9 . An inversion of the energy 
trend between these two structures, i.e. E^ e < can be seen for optB86b-vdW, PBE-D2 and 

PBE-TS, but not for PBE-D2, thus showing that the use of different van der Waals interactions can 
also affect results involving chemisorption. 

The inclusion of dispersion interactions in the optimisation of all the structures treated in this 
work allowed for the identification of the tilted structure (‘Tilted’ in Fig [I]). This was found to be 
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Figure 2 Electron density difference plots for the ‘Tilted’ (top row) and ‘Flat’ (bottom row) BTAH/Cu 
structures using three representative vdW-inclusive functionals. Green represents region of charge density 
depletion and red regions of accumulation. The isosurface level is 0.002 e/ajj for all structures. The distance 
between the azole N2 and the surface and between the centre of the benzene-like ring and the surface are given 
in A. 


the most stable configuration using optB86b-vdW (E ac j s = —1.07 eV) and degenerate to the ‘Flat’ 
structure with PBE-D2 (E a d s = —1.15 eV) and PBE-TS (E a d s = —1.29 eV). It is however unstable 
when optimised with PBE (the ‘Flat’ configuration is recovered) and unfavourable with vdW-DF, 
the adsorption energy being ~ 0.1 eV/mol higher than for the ‘Flat’ configuration. This result is 
in good agreement with previous work on functionalized benzene rings with PBE and vdW-inclusive 
functionals: it was seen that the inclusion of dispersion interactions is the driving element behind the 
chemisorption of molecules which would otherwise weakly physisorb using PBE 39 40 . From Fig. [ 2 ] it 
can be clearly seen that there is substantial rearrangement of the charge density at the azole end of 
the molecule and the Cu surface atoms to which it is bonded. This is indicative of the formation of 
a weak chemical bond between adsorbate and substrate. Therefore, from the results obtained for the 
‘Flat’ and ‘Tilted’ configuration with PBE and the vdW-inclusive functionals it is clear that dispersion 
stabilises these two configurations by bringing the molecule closer to the surface. This reinforces the 
chemical interaction of the molecule to the surface via the azole group, in the case of the ‘Tilted’ 
structure and via n bonding of the benzene ring, in the case of the ‘Flat’ structure. 

It can be noticed from the graph in Fig. [l] that the difference between the absolute adsorption 
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energies obtained with the functional yielding the weakest adsorption energies (PBE) and the strongest 
(PBE-TS) is ~ 0.9 eV/mol. In general, absolute adsorption energies will inevitably vary between 
exchange-correlation functionals^. However, studies of benzene on Cu(lll) have shown that the 
optB86b-vdW functional yields a binding energy in very good agreement with experiment, with the 
PBE-D2 and PBE-TS functionals overbinding and vdW-DF underbindingpl This makes optB86b- 
vdW the best candidate to study the relatively similar BTAH, and it has thus been chosen to obtain 
the results presented in the next sections. 

3.2 Coverage dependence of BTAH on Cu(lll) 
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Figure 3 a) Adsorption energy of BTAH on Cu(lll) as a function of coverage. Each point corresponds to a 
fully optimised adsorption structure. The filled symbols indicate a system where hydrogen bonding is present. 
The three coloured boxes define the three observed regimes: light blue for the physisorbed regime at low 
coverage, orange for the mixed adsorption regime at intermediate coverage and green for the chemisorption 
regime at the highest coverage, b-d) Top and side view of the most stable BTAH/Cu(lll) structure for each 
regime. In the top view the red box indicates the unit cell and periodic images are shown. In the side view the 
periodic images are not shown for clarity. Selected distances are given in A. b) Physisorbed regime: BTAH 
forms HB chains of flat molecules on Cu(lll). c) 1/8 ML coverage: mixed arrangement of HB flat and upright 
structures, d) 1/4 ML coverage: representative low-energy structure with vertically adsorbed molecules at a 
60° angle with respect to one another. 


Having established the most stable configurations for the isolated BTAH molecule on Cu(lll), the 

























stability of an extensive set of structures was evaluated with the optB86b-vdW functional as a function 
of coverage. The results, for 1/16 ML, 1/10 ML, 1/8 ML and 1/4 ML coverage, are shown in Fig. [3j 
Each point of Fig. [3^ corresponds to a fully optimized adsorption structure for BTAH on Cu(lll). 
Three regimes of adsorption can be observed: at low coverage (up to 1/10 ML), hydrogen-bonded 
physisorbed structures were found to be dominant, at intermediate coverage (1/8 ML) an adsorption 
mode consisting of a mixture of physisorption and chemisorption is found, and at high coverage (1/4 
ML) pure chemisorption is favoured. Hydrogen bonding is present in all low-energy configurations. It 
is dominant at low-coverage and it competes with other forms of intermolecular interaction at high 
coverage. 

In the range of pure physisorption, BTAH molecules form regular patterns in the form of hydrogen 
bonded chains with the molecules parallel to the surface with the N2 atom on the top copper site 
(Fig-Ht*)- The ‘Tilted’ structure in Fig.[l]has an adsorption energy ^ 0.40 eV/mol higher (less stable) 
than this HB configuration ( E ac j s = —1.46 eV/mol). This shows that the study of the interaction 
of isolated molecules on the surface gives only a limited description of the behaviour of the system 
and that the study of supramolecular aggregates is necessary to establish the full picture, even at low 
coverages. The adsorption energy of the HB chain at 1/16 ML can be decomposed into a contribution 
of ~ —0.92 eV/mol from the physisorption of BTAH onto the copper surface, and a ~ —0.59 eV/mol 
contribution from hydrogen bonding (Fig [5]). It is interesting to note that even with PBE, where the 
adsorption energy of a flat isolated molecule is very low, the flat hydrogen bonded structure is the most 
stable adsorption structure (with adsorption energy E ac js = —0.61 which is ^ 0.20 eV lower than for 
the ‘Up’ structure). Fig. [5] shows that this is due to the strength of the hydrogen bond (Ehb ^ —0.51 
eV/mol) which contributes to more than 90% of the adsorption energy. 

Hydrogen bonding dominates adsorption also at 1/8 ML, however the closer packing of the molecules 
makes steric repulsive interactions between the hydrogen atoms of the benzene-like rings relevant, and 
forces the system in a mixed flat/upright configuration (Fig. [ 3 J 3 ) . In these non-flat configurations, the 
highly-directional HBs become weaker (Ehb ~ 0.3 — 0.5 eV/mol according to the torsion) and other 
forms of interaction, such as N-Cu or n-n bonding, become relevant. The flat HB structure seen for 
1/16 ML and 1/10 ML, is still among the low energy structures (Fig. J3^t), with the chains sliding 
on top of one another because of the space restrictions. The 7r-bonds between the vertically stacked 
benzene-like rings add ^ —0.10 eV/mol to the total energy of the structure, however the interaction 
of the molecules with the substrate diminishes and the BTAH-Cu contribution reduces to ~ —0.67 
eV/mol on average. This results in an overall destabilisation of the HB chains, from E^ s = —1.46 
eV/mol at 1/16 ML to E a ^ s = —1.33 eV/mol. 

The same packing issues are relevant at 1/4 ML where the molecules adsorb in an upright configura¬ 
tion (a representative structure is shown in Fig. [3]d) . At this coverage, HB and non-HB configurations 
(in the form of HB dimers and upright molecules, respectively) coexist in the low-energy regime. In 
the HB structure (Fig. [4^), two upright BTAH molecules form a strong HB to the detriment of their 
interaction with the surface, with most N-Cu bonds much larger than the optimal length (^ 2. 00 A). 
In the non-HB case (Fig. §>), the BTAHs are preferentially packed at a 60° angle with respect to 
one another, in a trade-off between the most stable anchoring of the upright molecules to the surface 
(with N2 and H10 on a copper top site) and avoiding lateral steric interactions between the H atoms 
in the benzene-like ring. It can be seen in Fig. [3^i that, for 1/4 ML, there are two HB structure with 
E a ds ^ —1.00 eV/mol and they both correspond to HB chains (with one HB per molecule). The flat 
lying chain has already been discussed, while the chain of upright molecules owes its reduced stability 
to the distortion of the HBs whose strength for this configurations is only ^0.1 eV/mol. 

Fig. [5] highlights the differences obtained when optimising the low- and high-coverage structures 
with optB86b-vdW or PBE. It can be seen that high-coverage packing is less favourable with PBE than 
with optB86b-vdW. The difference in energy between the low- and high-coverage structures is ~ 0.3 
eV/mol with PBE, meaning that high-coverage configurations are higher in energy by almost 50%. On 


9 



Figure 4 Top (left-hand-side) and side (right-hand-side) views of the most stable structures for BTAH at 1/4 
ML coverage a) HB dimer structure, b) Non-HB upright structure. The red box shows the unit cell. Selected 
distances are shown in A. In the top views the unit cell and periodic replicas are shown. In the side view only 
the molecules in the unit cell are shown for clarity. 


the other hand, the difference with optB86b-vdW is ~ 0.2 eV/mol, less than 15% of the low-coverage 
interaction energy. Indeed, at high-coverage, where HB becomes less important, intramolecular in¬ 
teractions (vdW and electrostatic), shown in yellow in Fig. [5j become very important in stabilising 
the system. In the case of PBE the contribution of vdW is not accounted for and the high coverage 
systems are thus less stable. 

Overall we see that BTAH exhibits an incredibly rich phase behaviour when adsorbed on Cu(lll), 
with hydrogen bonding dominating low energy adsorption before packing effects kick in and non- 
H-bonded configurations become competitive. Moreover, significant differences between PBE and 
optB86b-vdW in the description of the adsorption behaviour of the molecules were seen for isolated 
molecules at low coverage and for high-coverage configurations. Van der Waals forces increase the in¬ 
teraction between the molecules and the substrate and between the molecules in the overlayer, meaning 
that more strongly bound and higher coverage overlayers can be obtained. This suggests that the use of 
vdW-inclusive functionals is necessary when studying the coverage dependence of molecular overlayers 
where intermolecular n bonding is possible. 
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Figure 5 Interaction energies for low coverage chains, high coverage dimers and upright molecules. The 
contributions to the interaction energy are represented by different colours: in blue are the contributions from 
the interaction of the molecules with the surface, in red interactions from hydrogen bonds and in yellow 
intramolecular interactions not from HB (van der Waals and electrostatic interactions). The components of the 
decompositions were calculated using Eqs. |5][6j 


4 BTA on Cu(lll) 


As already mentioned, the study of protonated BTAH is likely to be relevant for acidic environ¬ 
ments, whereas experimental evidence suggests that BTAH is deprotonated in alkaline solutions and 
in uhv hh6 H8] It 

has been suggested that BTA forms organometallic structures on copper surfaces, 
i.e. ordered networks of BTA molecules connected to one another via copper adatoms, in lieu of the 
lost hydrogen atom which allows the formation of hydrogen bonded networks (see e.g. 1118 and other 
references in 4 ^). Thus, understanding the structure of deprotonated benzotriazole on the Cu(lll) sur¬ 
face is important from the point of view of applications (in alkaline environments) and as a direct link 
and validation with experimental work performed in vacuum conditions. 


The most stable structures formed by BTA on Cu(lll) are presented in Sec. 4.1 and energetic 
considerations regarding the formation of the system, with respect to systems of fully protonated 

All calculations were performed with the optB86b-vdW 


BTAH molecules are detailed in Sec. |4.2 
functional. 
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Inhibitor 

77BTA 

^ads 

BTA (a,b) 

-3.36 

BTA-Cu ad (c) 

-3.94 

BTA-Cu ad -BTA (d) 

-4.02 

BTA-Cu ad -BTA stacked 

-4.24 

[BTA - Cu ad ] n (e) 

-4.55 

Corrosive atoms 

-®ads 

s 

- 4.79 

Cl 

-3.54 


Table 2 E^ A for the BTA/Cu(lll) structures shown in Fig. § All values were calculated using Eq. [I] and 
are in eV/mol. Also reported are adsorption energies of two typical corrosive atoms, Cl and S. These 
adsorption energies (also in eV) are with respect to gas phase Cl and S atoms. 

4.1 Structures and adsorption energies 

The adsorption of single deprotonated molecules on Cu(lll) and the formation of BTA-Cu a d organometal- 
lic complexes on the surface were considered. The most stable low-coverage (1/20-1/16 ML) structures 
are shown in Fig. [6j for a single BTA (panels a and b), for a BTA-Cu a d complex (panel c), for a 
BTA-Cu a d-BTA dimer (panel d) and for a [BTA-Cu a d] n chain (panel e). Their adsorption energies 
(calculated using Eq. [l] with BTA as a reference) are shown in the top part of Table [2j It can be 
seen that BTA interacts significantly more strongly (> 2.0 eV/molecule) with the surface than BTAH. 
No appreciable charge transfer was seen between the H atoms and the BTA molecule when mod¬ 
elled in the same unit cell, suggesting the molecule undergoes a dehydrogenation process (rather than 
deprotonation). 

For the isolated BTA, two degenerate low energy structures (E a d s = —3.36 eV/mol) were found: 
a heavily tilted structure (Fig. |6^), and an upright structure (Fig. [6 )d) . In the tilted structure the 
bonding to the surface is obtained via all three nitrogen atoms (resulting in a 0.05 A stretch in the 
N2N3 bond) and physisorption through the benzene-like ring. Only the N2-Cu bond has the optimal 
length of 1.96 A; N2 and N3 form weaker bonds with a length ^ 2.20 A. The upright structure is 
similar to the ‘Side’ structure for BTAH, and the interaction with the surface occurs via two N-Cu 
bonds. 

In the BTA-Cu a d complex (Fig. [6]^, E a d s = —3.94 eV/mol) the BTA molecule is tilted ~ 70° with 
respect to the surface and it forms three N-Cu bonds, two with the surface and one with the adatom 
(in the bridge site). Upright BTA-Cu a d structures were found to have higher energy by ^ 0.2 eV, 
although they were previously shown to be the most stable structure with pbeE) 

Thus, as with BTAH, we find that vdW forces alter the adsorption structures formed. Isolated and 
stacked (with a distance of ^ 4.4 A between successive copper adatoms) dimers were investigated. 

A number of almost degenerate low energy configurations (E a d s = —4.02 eV/mol) were seen for the 
isolated dimer, in contrast with PBE results where upright dimers with Cu ac j at the top side was found 
to be the most stable structure. The inclinations of the molecule with respect to the surface in these 
structures was found to range from upright to ^ 45° inclination, and the copper adatom was found 
to sit in either the top or bridge position (a representative configuration is shown in Fig. [6]T) . In 
general, the lowest energy configurations were found when the distance between the nitrogen atom in 
the BTA molecule and copper adatom is not larger than d(N-Cu a d) ^ 1.9 A and the distance of Cu a d 
with the substrate is not larger than d(Cu a d-Cu(lll)) ^ 2.5 A. In the case of stacked dimers, upright 
configurations were found to be preferred and the n — tv bonding between successive dimers stabilises 
the system by ^ 0.2 eV/mol. 

The strongest adsorption energy, E^ S A = —4.55 eV/mol, was found for the [BTA-Cu a d]n organometal- 
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lie chains. Fig. [6^ shows the most stable configuration, with the copper adatoms in hollow sites of 
the Cu(lll) surface and one vertical BTA, forming three N-Cu bonds, two with the adatoms and one 
with the substrate, and a flat lying BTA, forming two N-Cu bonds with the two copper adatoms. This 
structure is in good agreement with the structures found with PBE 8 10 with however a higher rotation 
of the molecule which, in our case, lies almost flat while with PBE it is tilted by ^ 60° with respect 
to the surface. Tests of stacked organometallic chains show that they do not significantly benefit from 
increased packing, yielding an adsorption energy which is degenerate with respect to the isolated chain. 
Indeed, the flat-lying molecule in Fig. [6^ has to rotate at a ^ 45° angle with respect to the surface in 
order to allow for stacking, and the tv bonding with the surface is replaced with n bonding with the 
benzene-like ring of the BTA in the parallel chains. 

When the adsorption energies for the BTA systems are compared to the adsorption energy of atoms 
like sulphur or chlorine (bottom of Table [2]), well-known active agents in the corrosion process of cop¬ 
per, it can be seen that their adsorption energies are very close, especially for chlorine. Although 
this is a fairly crude comparison (corrosive agents in a realistic system are likely to be sulphur- or 
chlorine-containing molecules rather than isolated atoms), these values provide a hint that competitive 
adsorption could be an important factor in the corrosion inhibition process of BTA on copper. 

4.2 Energies of formation 

A direct comparison between the stability of the systems discussed in this and in Sec. [3] is given by 
their formation energies, calculated using Eqs. [2]{3] for BTA and Ef orm =E a d s for BTAH (see Sec. [2]). 
Four scenarios for the formation of the BTA systems are shown in Fig. [T| (each point corresponds to a 
fully optimised BTA/Cu(lll) structure). In Fig. [7^ formation energies according to Eq. [2] are shown, 
i.e. it is assumed that the surface contains defects and that copper adatoms are already present on the 
surface (in the hollow site) before the adsorption of the BTAH molecules. The organometallic chains 
are the most stable system in this case. In Fig. [ 7)3 formation energies according to Eq. [3] are shown, 
i.e. it is assumed that the surface is initially atomically flat and that the adatoms are extracted from 
the bulk. The number of adatoms used in Eq. [3] is m = 0 for the BTA/Cu system, m — 1 for the BTA- 
Cu a d complex and for the BTA-Cu ac j-BTA dimers and m — 2 for the chains. In this case the stacked 
BTA-Cu a d-BTA structure have the lowest energy, ~ 0.2 eV/mol more stable than the organometallic 
necklaces. The two y axes report energies for the two systems formed by the dissociated hydrogen 
atom: adsorbed on the copper surface (left-hand-side) and forming H 2 molecules (right-hand-side). 

In Table [3] the formation energies for the most stable structures in all four cases are reported. 
These are compared to the formation energy of the most stable BTAH/Cu(lll) structure, the HB 
chains, and the blue highlighting indicates the deprotonated structures which are more stable than 
the BTAH/Cu(lll) system (Ef orm <E^^ H = —1.46 eV). When the formation energy of the copper 
adatom is not taken into account (top half of Table [5]) most deprotonated structures, except for the 
isolated BTA molecule, are found to be more stable than the protonated HB chains. If the formation 
energy of the copper adatom is considered (bottom half of Table [3]), only the organometallic chains and 
the stacked dimers are found to be more stable that the BTAH HB chains, and only if the hydrogen 
atom is assumed to adsorb on the surface. If the H atoms are assumed to associate into H 2 molecules 
the only competitive structure to the BTAH HB chains are the stacked dimers (Ef orm = —1.38 eV/mol). 
This shows that benzotriazole has a higher drive to deprotonate on defective surfaces with Cu adatoms 
than on atomically flat surfaces. 

The results obtained from these four extreme cases can be used to give an approximate description 
of the behaviour of a real experimental system. In the work of Grillo et al . 18 the copper surface is 
reconstructed and therefore the mobility of the surface copper atoms is low thus making this system 
close to the case in Fig. [ 7 ) 3 , i.e. the copper adatoms need to be extracted from the bulk. In this 
case (Table [3]) the stacked BTA-Cu a d-BTA dimers, with the dissociated H atoms adsorbed on the 
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Figure 6 Top and side view of the most stable deprotonated BTA structures on Cu(lll) at low coverage. The 
Cu adatom is shown in grey, a) Tilted configuration of BTA on Cu(lll). b) Upright configuration of BTA on 
Cu(lll). c) BTA-Cu a d complex. The most stable position for the copper adatom was found to be the bridge 
site, d) BTA-Cu a( }-BTA dimer. This is a representative structure of a low-energy configuration dimer with the 
copper adatom on the copper top site and a heavily tilted BTA. Other low-energy configurations include the 
copper adatom on the bridge site and upright BTA molecules, e) [BTA-Cu a d] n chain. The copper adatoms are 
found on the hollow sites. Selected distances are given in A and, for clarity, periodic replicas are not shown. 
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Structure 

r^H/Cu(lll) 771H2 

“*iorm ^form 

BTA (Fig. |6pb7 

BTA-Cu ad (Fig.g) 
BTA-Cu ad -BTA (Fig.g) 
BTA-Cu a d-BTA stacked 
[BTA — Cu ad ] n (Fig. §3) 

- 1.12 - 0.93 

-1.72 -1.51 

-1.82 -1.61 

-2.02 -1.81 

-2.31 -2.10 

Structure 

p H/Cu(lll),Cu ad p H 2 ,Cu ad 

"rorm "iorm 

BTA (Fig. g,b) 

BTA-Cu ad (Fig. gO 
BTA-Cu ad -BTA (Fig. 
BTA-Cu a d-BTA stack 
[BTA - Cu ad ] n (Fig. [( 

§i) 

ed 

IL 

-1.12 -0.93 

-0.86 -0.65 

-1.38 -1.18 

-1.59 -1.38 

-1.45 -1.24 


Table 3 inform for the BTA/Cu(lll) structures shown in Fig. [6] The top part corresponds to the data shown 
in Fig. UK The values are calculated using Eq. [2]z.e. the formation energy of the adatoms is not taken into 
account. The bottom part corresponds to the data shown in Fig. [TJd i.e. it is assumed that the copper adatoms 
are extracted from the bulk (Eq. [3|. The left column shows the formation energies when assuming that the 
dissociated H atoms adsorb on Cu(lll), the right column when assuming they form gas phase H 2 molecules. 
All values are in eV/mol. The values highlighted in blue are those where Ef orm < —1.46 eV, the formation 
energy of the BTAH flat HB structure. 


surface, are the lowest energy configuration with Ef orm = —1.59. This is in excellent agreement with 
the STM experimental results, where stacked dimers are observed at low coverage and where the 
reconstruction of the Cu(lll) surface, which is generally known not to reconstruct^, indicates the 
presence of impurities 42 (such as H atoms) on the surface. 

5 Discussion and conclusion 

Results for the structures and energies of BTAH and BTA on Cu(lll) have been presented here. 
The adsorption of fully protonated BTAH is relevant to applications where the environment is acidic, 
whereas the adsorption of BTA on copper is relevant to alkaline environments and to compare with 
experiments of BTAH deposition on copper in vacuum conditions, where the molecule is found to 
deprotonate. Benzotriazole is a complex molecule which requires treatment with a theoretical method¬ 
ology capable of simultaneously describing chemisorption and physisorption. In the present work DFT 
with a vdW-inclusive functional has been employed to optimise a large number of protonated and 
deprotonated structures on Cu(lll). 

We find that dispersion forces significantly alter the relative stabilities of adsorbed BTAH structures. 
In addition to this, in the lowest energy tilted structure an interesting interplay with chemical bonding 
is found wherein dispersion forces bring the molecule close to the surface thus enhancing the chemical 
bonding of the molecule with the surface via the triazole group. Whilst this is interesting and a 
potentially general effect, we have also shown that a single absorbed molecule provides limited insight 
into the behaviour of the overall system. Indeed, while isolated BTAH preferentially adsorbs on copper 
via the azole nitrogen atoms and 7r-bonding of the carbon ring with the surface, the lowest energy 
structure at low coverage are HB chains of flat-lying molecules. 

Overall we have found an incredibly rich coverage-dependent phase diagram for BTAH on Cu(lll). 
Three regimes were identified for the adsorption of BTAH on Cu(lll), as a function of coverage: a 
low coverage hydrogen-bonded regime, where the molecules preferentially adsorb flat on the surface, 
an intermediate regime, where mixed flat and upright structures are observed, and a high coverage 
regime, where the molecules adsorb upright. Steric interaction drive the change in configuration from 
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the flat-lying physisorbed to the upright chemisorbed configuration of the BTAH. 

The lowest energy configurations seen for BTA are either stacked BTA-Cu a d-BTA dimers or organometal- 
lic chains, according to whether the formation energy of the copper adatom is taken into account 
when calculating the stability of the complexes. Good agreement is seen with experimental results 
in UHV, where the molecule was found to adsorb via the azole moiety in a vertical or near-vertical 
manner 15 16 18 . In particular, the stacked dimer configuration was observed using STM to form on a 
reconstructed Cu(lll) surface, where the mobility of the copper atoms is likely to be low and therefore 
copper adatoms require a large amount of energy to form. In this case, the conditions of Fig. apply, 
and indeed the expected configuration from calculations are stacked BTA-Cu a d-BTA dimers. 

The use of a suitable exchange-correlation functional was found to be important for all the systems 
considered here. Indeed, for an isolated BTAH on Cu(lll) the ‘Flat’ configuration was found to be 
favourable only when vdW dispersion forces were accounted for, and not in the case of PBE. Moreover, 
the ‘Tilted’ low coverage structure, which is the most stable with the optB86b-vdW functional and is 
favourable with all vdW functionals tested, is instead unstable when optimized with PBE (the ‘Flat’ 
configuration is retrieved instead). At high coverage, the lack of any description of n-n bonding in 
PBE leads to a larger equilibrium distance between two BTAH molecules and thus to (comparatively) 
weaker adsorption energies for the high-coverage structures. In the BTA/Cu systems it has been 
seen that PBE favours upright adsorption, because of the lack of dispersion interactions between the 
benzene-like ring and the surface. When vdW interactions are accounted for a more complex behaviour 
is uncovered, especially for the BTA-Cu a d-BTA dimers where many degenerate low-energy structures 
are seen. 

Both BTAH and BTA offer the possibility of forming fairly close packed layers on Cu(lll) with 
little cost to the adsorption energy, and, for the case of BTA, with a gain in energy when the dimers 
are stacked, thus in principle offering a physical barrier to incoming corrosive molecules or atoms. 
There is however a large difference (~ 2 eV/mol) in the adsorption energy of the molecule with 
the surface between BTAH and BTA. Since benzotriazole performs better as an inhibitor in alkaline 
conditions, where the likelihood of deprotonation is higher, there might be a link between the strongest 
interaction of the molecule with the surface and inhibition. Indeed, the adsorption energy of BTA 
(~ 4.00 eV/atom) is fairly close to the adsorption energy of e.g. two well known corrosive agents for 
copper. It is likely that the actual corrosive agents in a corrosive solution are sulphur- or chlorine- 
containing molecules, rather than isolated atoms. However, our calculations on chlorine and sulphur 
atoms give a ballpark estimate which suggests that competitive adsorption could be the key here for 
the success of benzotriazole as a corrosion inhibitor. 
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Figure 7 Formation energies for an isolated BTA, for a BTA-Cu a d complex, for an isolated BTA-Cu a d-BTA 
dimer, for stacked dimers and for a [BTA-Cu a d] n necklace. The formation energies are calculated assuming 
that the H atom adsorbs on the surface (left-hand-side y- axis) or forms a gas-phase H 2 molecule 
(right-hand-side y- axis), a) Formation energies for a BTAH dissociating and adsorbing on a surface where 
adatoms are already presents, following Eq. [ 2 ] b) Formation energies for a BTAH dissociating and adsorbing 
on a surface where adatoms are extracted from the bulk, following Eq. [3] with m = 0 for BTA, m = 1 for 
BTA-Cu a d, dimer and stacked dimer, m — 2 for [BTA-Cu a d] n )- 
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